Exact General Relativistic Discs and the Advance of Perihelion 
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The advance of perihelion for geodesic motion on the galactic plane of some exact general rela- 
tivistic disc solutions is calculated. Approximate analytical and numerical results are presented for 
the static Chazy-Curzon and the Schwarzschild discs in Weyl coordinates, the Schwarzschid disc 
in isotropic coordinantes and the stationary Kerr disc in the Weyl-Lewis-Papapetrou metrics. It 
is found that for these disc models the advance of perihelion may be an increasing or decreasing 
function of the orbital excentricity. The precession due to Newtonian gravity for these disc models 
is also calculated. 
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The explanation of the anomalous precession of Mercury's orbit by Einstein in 1915 was one of the first successful 
predictions of General Relativity. In recent times the subject has received considerable attention with the possibility of 
high precision measurements of general relativistic effects in the orbits of binary pulsars like the PSR 1913+16 system 
discovered in 1974 [T3|. Schafer & Darmour [z| discuss in depth higher order general relativistic contributions to the 
C\| ' periastron advance of binary pulsars. Schafer & Wex [23[ calculated the periastron advance for a system composed 
of a Kerr black hole and an orbiting star. 

In the context of axially symmetric solutions of Einstein vacuum equations, Boisseau & Letelier Q calculated the 
O |. effect of different general relativistic multipole expansions in the advance of perihelion of test particles orbiting static 
axially symmetric attraction centres. Bini, De Paolis, Geralico, Ingrosso & Nucita [f| derived approximate expressions 
for the periastron shift for motion in static and stationary axially symmetric spacetimes. However, we did not find in 
the literature similar works for general relativistic solutions with matter, in p articular, discdike configurations. Several 
such solutions can be found in the literature, e.g., Mor gan & Morgan [Tjl; Bardeen & Wagoner Lynden-Bell & 
Pineault [13]; Chamorro, Gregory & Stewart Si; Lemos |l5j; Bicak & Ledvinka Q; Bicak, Lynden-Bell & Pichon j4j; 
Lemos & Letelier Neugebauer & Meinel [20l ]; Pichon & Lynden-Bell [2l[; Gonzalez & Espitia (lpj . 

The aim of this work is to study the advance of perihelion for motion of test particles in the galactic plane for a few 
exact solutions of Einstein field equations that represent disc-like configurations of matter (Bicak et al. 3; Gonzalez 
ff") \ & Letelier [ll|, [12]; Vogt & Letelier [U, HE 26]). We derive approximate expressions and also present numerical 
t-H ■ results. We find that when matter is present the periastron shift may be an increasing or an decreasing function of the 
orbital excentricity. The paper is divided as follows: in Section[n]we present the formalism to calculate the periastron 
• shift for relativistic elliptic orbits of test particles. This formalism is then applied in Section IlTTI to two exact models 
of static relativistic discs in canonical Weyl coordinates and to one solution expressed in isotropic coordinates. In 
Section HVl we calculate the periastron shift for a solution of a rotating disc obtained from the Kerr metric. In Section 
fVl we calculate the contribution of Newtonian gravity to the precession in the presented disc models so that it can be 
, separated from relativistic effects. Finally, in Section IVT1 we present a short discussion of the results. Along the work 
we take units such that c = G = 1 . 



II. ADVANCE OF THE PERIHELION IN RELATIVISTIC ORBITS 

In this section we derive the formulae to calculate the advance of perihelion for geodesic elliptic-like orbits in an 
axial symmetric spacetime with cylindrical coordinates (t, r, z, <p). We follow closely Bini et al. [5j. Let us assume a 
test particle is bound in an elliptic orbit on the plane z = 0. This orbit can be parametrized as 



1 + e cos x 

where d and e are, respectively, the ellipse's semi-major axis and excentricity, and x is a variable called relativistic 
anomaly. From equation JT]) we see that the minimum value r m = d(l — e) is obtained for x = and the maximum 
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value r p = d(l + e) when x = tt. At these points the equation dr/dy? that describes the shape of the orbit vanishes. 
The relation between the functions ip e x can be expressed as 



dip ed(l — e 2 ) sinx dip 
dx (1 + e cosx) 2 dr 



(2) 



r=r(x) 



where equation (|TJ) was used. By symmetry, the change in the coordinate (p when x decreases from it to is the same 
that when x increses from to 7r; thus the total change in the coordinate ip in one revolution is 2{ip{i:) — ip(0)), where 



ip(n) - ip(0) 



T ed(l — e 2 ) sinx dip 



dx- 



r=r(x) 



/ (1 + ecosx) 2 dr 

In a closed ellipse ip would change by 27r per revolution, so the orbit precesses by an angle 

Atp = 2(^(tt) - ip(0)) - 2tt 



(3) 



(4) 



in one revolution. In general it is not possible to express the integral equation ([3]) in terms of elementary functions; 
we will evaluate it numerically and also derive approximate expressions. 

III. ADVANCE OF THE PERIHELION AND STATIC RELATIVISTIC DISCS 

We study first the precession of perihelion for orbits in static relativistic disc models in Weyl coordinates and 
isotropic coordinates. 



A. Weyl Coordinates 

The general metric for a static axially symmetric spacetime in Weyl's canonical coordinates {t, r, z, ip) is given by 

ds 2 = -e 2 *dt 2 + e- 2 * [e 2A (dr 2 + dz 2 ) + r 2 dip 2 ] , (5) 

where <J> and A are functions of r and z. The Einstein vacuum equations for this metric reduce to the Weyl equations 
(Weyl [H,[29|) 



®,rr + — + <S>,zz = 0, 
r 

\2 *.2\ 



(G) 

A r =r(^-^), A z = 2r$ r $ z . (7) 
We shall consider two solutions of equations ©-([7]): the Chazy-Curzon solution (Chazyd; CurzonQ) 

(8) 



e 2 * = e ~ 2m/R 



e 2A = e -m 2 r 2 /fl 4 



where R = \/r 2 + z 2 , and the Schwarzschild solution, expressed as (Weyl 2 



$ = - In 
2 



i?i + R 2 - 2m 



Ri + R 2 + 2m 



A = - In 
2 



(i?i + R 2 f - 4to 2 



4i?ii? 2 



(9) 



with R\ = r 2 + (m + z) 2 , R\ = r 2 + (-m + z) 2 . 

Let us briefly recall a procedure to generate disc-like distributions of matter given a vacuum solution of Einstein 
field equations. Mathematically, it consists in applying a transformation z — > h(z) + a on a given vacuum solution and 
then calculate the resulting energy- momentum tensor using Einstein's field equations. Thin discs can be obtained if 
we choose h = \z\. For instance, Bicak et al. [H constructed thin discs using the Curzon solution equation §8§ and 
the Schwarzschild solution equation ([9]). On the other hand, thick discs can be constructed starting with the same 
vacuum solutions and using a class of even polynomials for h(z); see Gonzalez & Letelier [l2j and Vogt & Letelier 25 1 
for details. Also, a transformation originally proposed by Miyamoto & Nagai [l|| with h(z) = \/ z 2 + b 2 was used by 
Vogt & Letelier [2|| to generate relativistic disc-like distributions of matter from the Schwarzschild vacuum solution 
in isotropic coordinates. For our analysis, the advance of perihelion is always calculated on the galactic plane z = 0, 
where all the above mentioned transformations reduce to a constant. Thus, the results apply equally to thin and to 
thick discs. Henceforth this constant will be denoted a. 
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For timelike orbits on the z = plane, the Lagrangean associated to metric equation ([5]) reads 

2£ = -1 = -e 2 *i 2 + e 2 ( A -*)f 2 + r 2 e- 2 V, (10) 

where dots indicate differentiation with respect to proper time. Due to the independence of £ from t e ip, the conserved 
energy E and angular momentum h per unit mass can be introduced 

E = e 2 *t, h = r 2 e- 2 V (11) 

Using equation (fTTj) . the expression for the shape of the orbit follows from equation (fTOf 

1/2 

(12) 



dr r 

d(/? e A 



r 2 e~ 2 * (£ 2 e- 2 * - l) 



fc 2 



For an elliptic orbit with excentricity e and semi-major axis d, the two constants of motion can be calculated by 
substituting r m = d(l — e) and r p = d(l + e) in dr/dip = and solving the system. We have 

2 -2* p _ 2 ~2* m 2 2 -2(* p +* m )/- -2* ro _ p -2* p \ 

2 ^ 'pL- 'in'' ,2 _ 'p'm° r ° ,/ „n 

~~ r 2 p -4* p _ r 2 P -4<S> m ' — r ,2 p -4* p _ r 2 P -4<S> m ' \ 6 > 

where $ m = $(r m ) and $ p = $(r p ). 

To estimate the advance of perihelion for orbits in the z = plane for both disc models, it is reasonable to suppose 
that m/d and a/d are small quantities and expand equation ([3]) with equations (| 12[) — (| 1 3[) in multivariate Taylor series. 
For the Curzon disc, the expansion up to third order reads 

6ttto 3?ra 2 7rm 2 (44-9e 2 ) 67rma 2 (6 + e 2 ) 7rm 3 (192 - 53e 2 ) 

d{l-e 2 ) ~ d 2 {l-e 2 ) 2 + 2d 2 (l-e 2 ) 2 d 3 (l - e 2 ) 3 + 2d 3 (l-e 2 ) 3 ' ( ' 

while for the Schwarzschild disc we obtain 

6?rm 37ra 2 37rm 2 (14 - 3e 2 ) 6nma 2 (6 + e 2 ) 37rm 3 (56 - 19e 2 ) 

^~ d{l-e 2 ) ~ d 2 {l-e 2 ) 2 + 2d 2 (l-e 2 ) 2 d 3 (l - e 2 ) 3 + 2d 3 (l - e 2 ) 3 ' (15) 

Both expansions with a = agree with those presented by Bini et al. [H up to second order. In equations (|14[) - ([T5|) 
the terms corresponding to the vacuum solutions are all positive, whereas the ones related to the presence of matter 
(parameter a) have negative signs. When matter is absent, the angle of advance is an increasing function with respect 
to excentricity, but expansions (|14p - (|15|) suggest that this may not be true in the present disc models. By imposing 
d(A(p)/de = in equations (fT4")) - (|15[) . the following expressions for the parameter a are obtained, respectively, 



.2 



a 



m 12d 2 (l - e 2 ) 2 + md{\ - e 2 )(79 - 9e 2 ) + m 2 (523 - 106e 2 



(16) 



4[d(l - e 2 ) + m(19 + 2e 2 )] 
These expressions, evaluated for e = and e = 1, yield the results 



12[d(l-e 2 )+m(19 + 2e 2 )] 

i 2 fU9-3Rft 2 Yl 

(17) 



2 _ to [4d 2 (l - e 2 ) 2 + md(l - e 2 )(25 - 3e 2 ) + m 2 (149 - 38e 2 )] 



, m(12d 2 + 79md + 523to 2 ) 2 139to 2 



a = 12(d+19m) ' a = IT (CurZ ° n) ' (18) 

9 m(4d 2 + 25md + 149to 2 ) 2 37to 2 ,„ . s , ^ 

a= WTmn) ' a =^T (Schw - } - (19) 

Equations ()18p ~ (|19[) give an estimate for the ranges in the parameter a for which the angle of advance as a function 
of excentricity has a critical point. 

Figs QJa)-(c) show the angle of precession as function of excentricity for the Chazy-Curzon disc with parameters 
m/d = 0.01, a/d — 0.05 in fig. QJa), a/d = 0.075 in (b) and a/d = 0.1 in (c). The curves with solid lines were 
calculated by numerical integration of equation ([3]), those with dashed lines represent expansion equation (fT4^) with 
terms up to second order, and the dotted lines the same expression but with terms of third order. Figs [DJa)— (c) show 
the results for the Schwarzschild disc with the same value of parameters. In both cases the numerical results and the 
expansion up to third order agree well for small values of e. Using the estimates given by equations (|18p and (jTT 
the intervals with critical points would be 0.0128 < a/d < 0.0948 and 0.0115 < a/d < 0.0946, respectively. 
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FIG. 1: The advance of the perihelion A<p as function of excentricity e for the Chazy-Curzon disc. Parameters: m/d = 0.01, 
a/d = 0.05 in (a), a/d = 0.075 in (b) and a/d = 0.1 in (c). Solid lines: numerical integration of equation (|3}. Dashed lines: 
values obtained from equation (|14| up to terms of second order. Dotted lines: the same expansion with terms of third order. 
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FIG. 2: The perihelion shift Aip as function of excentricity e for the Schwarzschild disc. Parameters: m/d — 0.01, a/d = 0.05 
in (a), a/d — 0.075 in (b) and a/d = 0.1 in (c). Solid lines: numerical integration of equation (J5J. Dashed lines: values obtained 
from equation (|15p up to terms of second order. Dotted lines: the same expansion with terms of third order. 
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B. Isotropic Coordinates 



The line element in isotropic form in cylindrical coordinates (t, r, z, if) may be expressed as 

ds 2 = -e 2 *di 2 + e 2A (dr 2 + dz 2 + r 2 d<^ 2 ), 

where the $ and A are only functions of r and z. The vacuum Schwarzschild solution for metric equation (|20[) has 
the form 



(20) 



2* 



1 2R 



2R 



■2A 



m \' 
2RJ 



(21) 



where R = \/r 2 + z 2 . Also here disc-like distributions of matter can be generated by app lying co nvenient transfor- 
mations on the z coordinate (see, for example Gonzalez & Letelier [T2I; Vogt & Letelier [24. [251 126|) . as was discussed 
in Section IrTLA"! 

For metric equation (|20)) . the shape of the orbit of a test-particle confined on the z — plane is described by 



dr 
dtp 



2 D 2A (£2 e -2* 



re 



h 2 



1/2 



The constants of motion E and h are now given by the following expressions 



E 



ry>2 p2Ap „2 p2Ap 



„2 2 2(A P +A m ) 
' p 1 m c 



-2$„ 



-2*„ 



r 2g2(Ap-* p ) _ r 2 e 2(A m -* m ) ' 



T ,2 p 2(Ap-$p) _ r 2 „2(A,„-$ m ) 
' p & ' m 



(22) 



(23) 



with the same notation as defined in Section IIII Al 

For the Schwarzschild disc in isotropic coordinates, an approximate expression for the precession of perihelion reads 



6nm 



37Tfl 2 



37TTO 2 (14-3e 2 ) 6nma 2 (6 + e 2 ) 3nm 3 {57 - 16e 2 ) 



d(l-e 2 ) d 2 (l-e 2 ) 2 



2d 2 (l-e 2 ) 2 



d 3 (l-e 2 ) 3 



2d 3 (l-e 2 f 



(24) 



Comparing equations (fT5|) and ([24| reveals that they are almost identical, the difference beginning only in the last 
term. The calculation of d(A(p)/de — provides 



in 



[4d 2 (l - e 2 ) 2 + md(l - e 2 )(25 - 3e 2 ) + m 2 (155 - 32e 2 



4 [d(l - e 2 ) + m(19 + 2e 2 )] 
which evaluated on e = and e — 1 gives, respectively 

2 m(Ad 2 + 25md+ 155m 2 ) 



4(d+ 19m) 



41m 2 
28 



(25) 



(26) 



Figs (Ha)-(c) display some curves of the angle Aip as function of the excentricity e for the Schwarzschild disc in 
isotropic coordinates with parameters m/d = 0.01, a/d — 0.05 in fig. [3]Ja), ajd — 0.075 in (b) and a/d = 0.1 in (c). 
The curves with solid lines were calculated by numerical integration of equation ((3|) , those with dashed lines represent 
expansion equation (|24|) with terms up to second order, and the dotted lines the same expression but with terms of 
third order. For this example, equation (j2"6")) estimate an interval 0.0121 < a/d < 0.0947. Finally in figs IUa)-(c) the 
three disc models are compared. Parameters were taken m/d = 0.15 e a/d = 0.175 in fig. [4ja), a/d — 0.225 in (b) 
and a/d = 0.275 in (c). All curves were obtained by numerical integration of equation ©. Solid lines represent the 
results for the Chazy-Curzon disc, dashed lines for the Schwarzschild disc in Weyl coordinates and dotted lines for 
the Schwarzschild disc in isotropic coordinates. Remember that as the parameter a is increased, all the discs become 
less relativistic (Bicak et al. 3 Vogt & Letelier I24T) . This is reflected in the numerical values of the precession angle, 
which are greater in fig. 0Ja) than in (c). Qualitatively the curves for the three disc models are similar. As suggested 
by the expansions equations (| 1 5(1 and (|24jl . both models obtained from the Schwarzschild solution give quite similar 
results for low excentric orbits. 
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FIG. 3: The advance of the perihelion Aip as function of excentricity e for the Schwarzschild disc in isotropic coordinates. 
Parameters: m/d — 0.01, a/d = 0.05 in (a), a/d = 0.075 in (b) and a/d = 0.1 in (c). Solid lines: numerical integration of 
equation ([3j|. Dashed lines: values obtained from equation (|24| up to terms of second order. Dotted lines: the same expansion 
with terms of third order. 




FIG. 4: The advance of the perihelion Aip as function of excentricity e for the three disc models. Parameters: m/d = 0.15, 
a/d — 0.175 in (a), a/d = 0.225 in (b) and a/d = 0.275 in (c). Solid lines: Chazy-Curzon disc. Dashed lines: Schwarzschild 
disc in Weyl coordinates. Dotted lines: Schwarzschild disc in isotropic coordinates. 
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IV. ADVANCE OF THE PERIHELION AND STATIONARY RELATIVISTIC DISCS 



In this section we investigate the effect of rotation on the perihelion shift for stationary relativistic discs, in particular 
discs generated from the vacuum Kerr metric. We begin with the metric for a stationary axially symmetric spacetime 



ds 2 



e 2$ (dt + Adip) 2 + e- 2 * [r 2 d^ 2 + e 2A (dr 2 + dz 2 )] , 



(27) 



where $, A and A are functions of r and z. The vacuum Kerr solution for metric equation (|27p may be written as 

(28) 
(29) 
(30) 



lln 

2 




R2) 2 


- 4m 2 + a 2 {Ri - 


R2?lo 2 \ 


. (Ri- 




2m) 2 + a 2 {R 1 - 


R2YI0 2 \ 


Im 

2 


"(i?i4 


R2) 2 


~ Am 2 + a 2 (R 1 - 


# 2 )> 2 






4i?ii? 2 




am 


+ R2- 


f 2m) [4ct 2 - {R x 


- R2) 2 ] 


a 2 


(i?i4 


R2) 2 


- Am 2 + a 2 {Ri - 


-R2?lo 2 ' 



where m and a are, respectively, the mass and the Kerr parameter, R\ = d r 2 + (z + a) 2 , Ri = \J r 2 + (z — a) 2 and 
u = \/m 2 — a 2 . Following the same procedure taken in Section HTT1 the orbit's shape of a test particle on the plane 
z = for metric eq. ([27]) is described by 



dr 

dip 



r 

7a 



1) 



1 1/2 



- 1 



(EA + hy 

The conserved energy E and angular momentum h are found by solving the system of equations 



2 -2* 
r p e 



r 2 e 



" (E~e 

e 



2-2$,, 



2<S> m ^£j2 „- .!•!• 



1) - (EA p + hf = 0, 

- 1) - (EAm + h) 2 = 0, 



(31) 



(32) 
(33) 



where A p = A(r p ), A m = A(r m ), and where again the same notation was used as defined in Section UlI Al For a given 
excentricity e and semi- major axis d the system of equations admit two distinct solutions, corresponding to prograde 
(h > 0) and retrograde (h < 0) orbits. After applying a convenient transformation on the Kerr solution e qua tions 
(|28|) -l[50 |) to generate stationary disc-like distributions of matter (see Gonzalez & Letelier fill: Vogt & Letelierl27t) . we 
assume the ratios m/d, a/d e a/d to be small and expand equation ([3]) in scries up to third order 30] 



6nm 



Z-ko 2 



37rm 2 (14-3e 2 ) 



3?ra 2 



127ram 3 / 2 (5 - e 2 ) 



d(l -e 2 ) ^ d 3 / 2 (l-e 2 ) 3 / 2 d 2 (l-e 2 ) 2 2d 2 (l-e 2 ) 2 rf 2 (l - e 2 ) 2 ^ d 5 / 2 (l - e 2 ) 5 / 2 

67rma 2 (6 + e 2 ) 37rm 3 (56 - 19e 2 ) 67rma 2 (12 - e 2 ) 



d 3 (l-e 2 ) 3 2d 3 (l-e 2 ) 3 d 3 (l-e 2 ) 3 

where the minus (plus) sign refers to prograde (retrograde) orbits. The solution of d(Aip)/de = yields 

1 r r . .0 2 \2 



(34) 



4[d(l-e 2 ) + m(19 + 2e 2 )] 



{m [4d 2 (l - e 2 ) 2 + md(l - e 2 )(25 - 3e 2 ) + m 2 (149 - 38e 2 ) + 4a 2 (35 - 2e 2 )] 



+Aad 1 / 2 {\ - e 2 ) 1 ' 2 [ad^ 2 (l - e 2 ) 1 ' 2 T 2m 1 ' 2 d{\ - e 2 ) T m 3/2 (23 - 3e 2 )] } , (35) 
which for e = and e = 1 simplifies to 



i(4d 2 + 25md + 149m 2 + 140a 2 ) + Aad 1 / 2 (ad 1 / 2 =p 2m 1 / 2 rf =p 23m 3 / 2 ) 

4(rf+ 19m) 

, 37m 2 + 44a 2 

a = , 

28 



(36) 
(37) 



respectively. 

Figs EJa)-(b) show the advance of the perihelion as function of excentricity for the Kerr disc with parameters 
m/d — 0.01, a/d — 0.075 and a/d = 0.001. The curves in (a) are the results for prograde orbits and (b) for 
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FIG. 5: The perihelion shift Ay> as function of excentricity e for the Kerr disc. Parameters: m/d — 0.01, a/d = 0.075 and 
a/d = 0.001. Fig. [5|a) displays results for prograde orbits and (b) for retrograde orbits. Solid lines: numerical integration of 
equation ((3| . Dashed lines: values obtained from equation (|34[) up to terms of second order. Dotted lines: the same expansion 
with terms of third order. 




FIG. 6: The perihelion shift Aip as function of excentricity e for the Kerr disc. Parameters: m/d — 0.15, a/d = 0.05, 
a/d — 0.175 in (a), a/d — 0.225 in (b) and a/d = 0.275 in (c). Solid lines: disc without rotation. Dashed lines: prograde orbits. 
Dotted lines: retrograde orbits. 



retrograde orbits. Curves with solid lines were calculated by numerical integration of equation (|3|), those with dashed 
lines represent expansion equation (|34[) with terms up to second order, and the dotted lines the same expression 
but with terms of third order. For these values we obtain from equation (|57|) intervals 0.0116 < a/d < 0.0970 and 
0.0116 < a/d < 0.0972 for prograde and retrograde orbits, respectively. The results of numerical integration of the 
exact expressions are depicted in figs EJa)-(c) with parameter values m/d = 0.15, a/d — 0.05, a/d = 0.175 in (a), 
a/d = 0.225 in (b) and a/d = 0.275 in (c). The values for prograde orbits are represented by dashed lines, retrograde 
orbits by dotted lines, and solid lines the perihelion shift without rotation (Weyl disc). We note that for prograde 
orbits the Kerr parameter lowers the angle of precession and has an opposite effect for retrograde orbits. The signs 
in the expansion equation (|34|) also predict these effects. 

Finally in figs[7{a)-(b) the parameters m/d = 0.15 and a/d — 0.225 were held constant and the Kerr parameter was 
changed: a/d = (solid lines), a/d = 0.05 (dashed lines) and a/d = 0.1 (dotted lines). Prograde orbits are shown in 
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(a) (b) 
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FIG. 7: The perihelion shift A<p as function of excentricity e for the Kerr disc. In (a) are shown prograde orbits and in (b) 
retrograde orbits. Parameters: m/d = 0.15, a/d = 0.225, a/d — (solid lines), a/d = 0.05 (dashed lines) and a/d — 0.1 
(dotted lines). 



fig- EI a ) and retrograde orbits in (b). Rotation has the same effect on both types of orbits as observed in fig. [SJ 



V. EFFECT OF NEWTONIAN PRECESSION 



The advance of perihelion calculated for the relativistic models of discs presented in Sections Mil and IIVI has two 
origins: one due to relativity and other from purely Newtonian gravity, since any flattened body will generate a 
perihelion precession. Thus, it would be interesting to separate the relativistic from the Newtonian contributions to 
the precession. 

It can be shown that in the non-relativistic limit the above mentioned disc models reduce on the z = plane to the 
Kuzmin model [l4| 

m 

Vt + a z 

The orbital equation in the usual cylindrical coordinates reads 

>\ 1 1/2 

(39) 



dr 



2r 2 (E M - $) _ 1 1/2 



h 2 

where Em is the conserved mechanical energy of the test particle 



4% - r^ m , , 2r 2 r 2 m (<5> p - <5> m ) 
E M = 2 _ 2 ' and h = — TZ~^2 > ( 40 ) 

^ p ^p 

whith the notation as defined in Section IIII Al Proceeding as in the previous Sections, it is possible to deduce 
approximate expressions for the perihelion shift. Using equations ([3]) and ([38|) - (|40|) , one has the following expansion 

^ = -^0^2+O{{a/dn (41) 

Note that the angle of advance is independent of m, which is canceled in the fraction in equation (|39| . Thus, the first 
term of equation (|4"Tj) , which also appears in expansions (Ti~4")) . (TT5)) , (|24[) and (f3"4"]) , is the purely Newtonian contribution 
up to third order to the perihelion shift. 

As a numerical example, in figs [SJa)— (c) we display curves of Aip as function of the excentricity e for the 
Schwarzschild disc in isotropic coordinates with m/d — 0.01, a/d — 0.175 in (a), a/d = 0.225 in (b) and a/d = 0.275 
in (c); the same values that were used in fig. 2J The solid lines represent the total angle of precession; the dotted 
lines represent the precession due to Newtonian effects obtained from the numerical integration of equation ([3]) with 
equations (|38|) ~([40 |) . and the dashed lines are the difference between the previous two values. It is seen that the pre- 
cession rate due to Newtonian gravity is in the opposite sense to the relativistic precession. Also for less relativistic 
discs the Newtonian contribution is more significant, as expected. 
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FIG. 8: The different contributions to the perihelion shift Atp as function of excentricity e for the Schwarzschild disc in isotropic 
coordinates. Parameters: m/d = 0.15, a/d = 0.175 in (a), a/d = 0.225 in (b) and a/d = 0.275 in (c). Solid lines: the total 
shift. Dotted lines: the purely Newtonian contribution. Dashed lines: the relativistic contribution. 



VI. DISCUSSION 



We studied the advance of perihelion for elliptic orbits of test particles in geodesic motion on the galactic plane 
for relativistic static and stationary disc models. We derived approximate expressions for the perihelion shift and 
compared them with the numerical integration of the exact solutions. The results show that the angle of advance 
can increase as well as decrease with increasing excentricity. We have that highly relativistic discs favours the first 
situation, and the advance of perihelion decreases with excentricity when the discs become less relativistic. The effect 
of rotation was also studied for a particular stationary disc model based on the Kerr solution. We found that the 
Kerr parameter lowers the perihelion shift for prograde orbits and increases it in the case of retrograde orbits. We 
also calculated the different contributions (Newtonian and relativistic) to the advance of perihelion for the relativistic 
disc models. 

Our conclusions are based on the study of a few exact general relativistic disc models, in particular, the Miyamoto- 
Nagai model, that presents some characteristics of real galaxies. We believe that the results found may be common 
to other models of galaxies. 
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